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Abstract 

The behavior of a non-spherical particle in a viscous, plane channel flow is studied by means of 
a combination of analytical technique and geometrical reasoning. An efficient implementation of 
Lamb’s general solution is adopted, allowing rapid calculation, in the perturbative quasi-spherical 
regime, of forces and torques coming from the particle interaction with the flow gradients and from 
the presence of the walls. Particular attention is given to the role of symmetries in the problem. 
A classification of the effects of departures from spherical shape on the motion of a particle in 
various kinds of shear flow is provided. Additional applications to problems involving particles 
with streaming boundaries are presented. 
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I. Introduction 


The motion of a spheroidal particle in a viscous shear flow is one of the classical problems of 
fluid mechanics, with important applications in diverse fields such as deposition and sedimentation 
processes , hemodynamics |§, , swimming of micro-organisms |Q, ||, D and the study of diffusion 

in suspensions It is well known that the simplest suspension model, in which the particles 

are rigid and spherical, and the flow is purely viscous, does not account for important physical 
processes such as the transverse migration of particles with respect to the flow lines. This is 
basically a consequence of the symmetry properties of the Stokes equation which describes the 
dynamics of a fluid in creeping flow conditions. 

Different physical mechanisms leading to the symmetry breaking necessary for transverse mi¬ 
gration have been analyzed in the literature j|]. Besides the effect of non-spherical shapes, the role 
of inertia was considered, both in the presence of walls bounding the fluid jl^, and in inflnite 
space situations (see 0 for recent references), while, Chan & Leal |Q, analyzed the role of 
non-Newtonian corrections to the constitutive law of the suspending fluid. 

The velocity disturbance produced around a non-spherical particle is very difficult to calculate, 
and, with a few exceptions such as the ellipsoid in a linear shear (see e.g. 1^), no analytical 
solutions exist. For shapes close to spherical, however, perturbation theory becomes possible and 
analysis of the modification to the particle dynamics produced by different kinds of deviation from 
sphericity and of external flow, can be carried on. 

An approach to the establishment of a perturbation theory is the one used e.g. in 0 and in 
1^, based on the use of the reciprocal theorem 00. One of the advantages of this approach 
consists of the fact that no information on the perturbation produced by the particle in the given 
fluid flow is necessary. Rather, the simpler problem of the translation (or rotation) of an identically 
shaped particle in a quiescent fluid must be considered. 

Direct perturbation expansion of Stokes equation, using Lamb’s general representation [0 , 
however, has the advantage of highlighting the symmetry structure of the fluid perturbation and 
its relation to that of the particle shape. The problem is that a grater amount of information seems 
to be necessary in order to determine the forces and torques necessary to obtain the particle law 
of motion. 

It turns out that it is possible to cast this representation into a form requiring an amount of 
perturbed velocity fields calculation, equivalent to that of the reciprocal theorem approach (see 
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also In the following sections, it will be shown how this can be done using an approach that 

is basically an adaptation of standard quantum mechanical matrix techniques. In this approach, 
the forces, torques and lift velocities can be obtained, in the spherical case, as components of the 
external flow velocity on the particle boundary and of possible streaming velocities of the particle 
surface. In the non-spherical case, calculating such forces and torques results in the computation of 
matrix elements involving these basis functions and the unperturbed velocity field (plus streaming), 
but not the velocity perturbation by the particle. 

In the next two sections Lamb’s representation will be recalled and the formalism for the 
calculation of forces and torques will be derived. In section IV, some applications to problems 
involving spherical particles with moving boundaries will be described; in particular the calculation 
of a bound for the swimming efficiency of micro-organisms, which is lower than the one obtained 
in iP, will be given. In section V, the behavior of a spheroidal particle in a plane channel flow 
will be considered, separating the various contributions from interaction with the walls and with 
the flow gradients. Various geometrical proofs on the functional form of forces and torques and 
on their relation to symmetry properties of the particles will be presented in the general strongly 
non-spherical case. Section VI will be devoted to the conclusions. Details of the calculations of 
forces and torques, and of the effect of the walls will be presented in the appendices. 


II. Revie-w of Lamb’s representation 


In a perturbative analysis around a spherically symmetric situation, it is convenient to expand 
the velocity field in terms of vector spherical harmonics ^ =s,l,m: 


'"W = + vf^{x)Yf^{xi) + vZ{x)Y^{xi)], (2.1) 

Im 

where Xi = Xi/x are the cosines of x in an appropriate reference frame and the superscripts {sem}, 
standing for scalar, electric and magnetic, come from the origin of this basis as a tool in the study 
of electromagnetic waves (see [ p^ ) In terms of standard spherical harmonics the functions are 
defined as 




YL(^.) 


x\7Yim{xt) 

\/l{l + 1) 


v™ ^ Y]Yim{Xi) 

rry: - 


( 2 . 2 ) 


where x = x ^x. It is natural to adopt a standard bra-ket quantum mechanical notation for the 
calculation of components and matrix elements: \^lm) = Y()^, {film\ = Y()0 and also: 


{ni 


TO V) = 


dfl; 




/■(x) and {film 


l/|v) = J I 


d^^^Y(;;)(xO • v(x)/(x) 


(2.3) 
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Here, dfia, is the solid angle differential around x and the left hand side still depends on the radial 
coordinate x. With the definition provided by Eqn. (2.2) (which differs from the one adopted 
in the vector spherical harmonics result normalized to one: = 1. Hence, one 

has the standard expression for velocity components v^^{x) = {^lm\v). These components can be 
expressed in terms of plain spherical harmonics integrals (see Appendix A): 

[= / ^^xYi*^{x^)± ■ v(i?x) 

\ vfmiR) = W + l))-i J dn,Y;^ix.)[3± ■ v(i?x) - V • (xv(Ax))],=« (2.4) 

Urm(^) = -W + 1))“^ /df2xh));,(i*)V • [v(i?x) X 

At stationarity, an incompressible fluid in creeping flow conditions, obeys the Stokes and continuity 

equations: 

= VP and V • V = 0, (2.5) 


where /i and P are the dynamic viscosity and the pressure. In terms of vector spherical harmonics. 
Lamb’s representation for the solution to eqns. (2.3-4) can be cast into the form: 

Vfm. = airaX~^ + blmX~'^~\ 

=aimX~^ - \[j^bimX~'^~\ (2.6a) 


'^Im 

''^tm — 


V^/(/+l) 


Sl‘^=ClmX ^ 


in an external domain, and: 


' + blmX’' \ 


,e _ 1+3 


VT^ = 


i/i(i+i) 




( 2 . 66 ) 


Xi^ = XlmY. 

in an internal domain. One finds also for the pressure field, respectively in an external and internal 
domain: 


V^ 4/-2 

-P = M , aimX 


Im 


i + i 


Yim and P — /r 2 ^ , x^imX Yi^n- 


Im 


I 


(2.7) 


and for the vorticity u) = V x v, in an external domain: 


-J2i~ VWd^ci^x ^ -h IcimX ^ - ^=^=aimx ; (2.8a) 

^ vW+T) ) 


iV = 

and in the internal domain: 
1 


UJ = 


-J2(- x/Z(7 + T)q„x'YL -{1 + l)cimx‘Y<i^ + -=^a,„x'+iY” ). (2.86) 

The velocity field is thus partitioned into contributions from the pressure (the aim terms), and into 
purely potential flow (the 6/m terms) and pressure free, vortical contributions (the c/m terms). 


4 














A sequence of formulae of practical use can be derived. Expressing the components of the 
velocity held in terms of velocity boundary conditions at a reference distance x = R: v(i?x) = 
U(x), one has, in an external domain: 


-L'i (2- 

.(I) 

while, in an internal domain: 


0 + (2 - 0 (f) ') UL + yWTT) (z - (f) ') t/fj (f)' 


(2.9a) 


^Im 


= 2 


m _ rrm 
^Im ^ Im 


(-(z + i) + (z + 3)(|) yf^ + yi(TTT)(i-(^§) )c/L](f) 

' (i+ 3) v^( -1+(f) +0+3 +1) (^)') t 

(t)' 


Z+1 


Z+1 


(2.96) 


Similarly, it is possible to express the force on a surface bounding the huid in terms of the velocity 
components on that surface. The force density f exerted by the huid on a spherical surface element 
at X = R, is given by the normal component of the stress tensor Sij = —PSij + atj: f = ±x • S, 
where (7^ = fi{diVj + djVi) is the viscous stress and the ± signs refer to the case of the surface 
bounding the huid, respectively, from the inside and from the outside. Using Eqns. (2.7) and 
(2.9a-b), one hnds, in the case of an external domain (see also plf): 


fs — M. 
J Im R 

fe _ ii 
J Im R 


21 +3Z+4 rrs i o / L tt' 
/+1 ^Im '^\l /-LI 


I -L 


l-\-l ^ Im 

Uf„-(2Z + l)C/L 


(2.10a) 


=-L(Z + 2)t/,- 


and 


{ fs _ A£ 

Jim R 

fe = ii 

Jim R 

in that of an internal domain. Inverting Eqns. (2.10b) and substituting into Eqns. (2.9a-b), it 
is possible to obtain expressions for the velocity held produced by a distribution of forces on the 
surface x = R. 


2R+1+3 TTS 
I ^Im 


. q, /hhir/' 

I '^Im 


- (2Z + 1)[/L 


(2.106) 


III Representation of forces and torques on particles 

The force F exerted on a particle by the huid is obtained integrating the expression for the force 
density provided by Eqns. (2.10a-b) on a spherical surface surrounding that particle. As it is well 
known ||r^ , the only term contributing to the force is the one coming from the aim components of 
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the velocity perturbation [see Eqn. (2.6a)]: 

Fi = 2 ^'/^Reaii; F 2 = F 3 = — 2 ^y/§Traio, (3-1) 

Now, the coefficients aim can be expressed in bra-ket notation as: an = |{(sll| +-\/2(ell|}|v) and 
similar for m ^ Thus, the force on the particle can be obtained calculating the component of the 
perturbed velocity on adequate basis functions, which are combinations of ’’scalar” and ’’electric” 
spherical harmonics, and which have the expression: 


|Fi) = ^a:[-3/2,0,0]; jF^) = ^^[O,-3/2, 0]; \Fs) = fix[0,0,-3/2]- 


(3.2) 


The torque T can be calculated in a similar way, starting from the expression R = f dfixSijkXjfkX^ 
= f dflxCijkXjXiaik and this leads to the result that now the ’’magnetic” I = 1 components of the 


velocity perturbation are to be taken in exam 117 : 

Ti = —T2 = dV^filmcn', T3 = 2 -\/^/rcio 


(3.3) 


In a way analogous to the case of the force, the torque can then be calculated extracting the 
appropriate components of the velocity field. The associated basis functions are proportional to 
’’magnetic” harmonics with ^ = 1: 


jTi) = 3fix[0, X 3 , -X 2 ]] IT 2 ) = 3 iix[x 3 , 0, -xij; jTs) = 3^x[x2, -xi, Oj. (3.4) 

In principle, Eqns. (3.2) and (3.4) require knowledge of the velocity perturbation around the 
particle to calculate force and torque components, however, for particle shapes close to spherical 
this is not the case. In the case of a radius R spherical particle, no-slip boundary conditions on 
the particle surface imply: 

U(x) = U(x) -f U(x) (3.5) 

where U and U are the value &t x = R oi the external flow v and of the velocity perturbation 
V, as seen in the reference frame translating with the particle, and, apart of the components 
representing a bulk rotation, U is non-zero only if there is streaming of the particle surface (or 
boundary conditions which allow for the presence of slip). Hence, the force and the torque on this 
particle will be given by the expressions: 

F = (F|{U-tj}) and T = (T|{U - U}) (3.6) 

where (F| and (T| are the vectors whose components are the ’’bras” (Fij and (F^j. (Clearly, in 
the absence of torque, = U^)- Thus, trivially, no information on the velocity perturbation 
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produced by the particle in the fluid is necessary to calculate forces and torques on the particle. 
This, independently of possible complication in the form of the surface streaming U, which makes 
this approach useful in problems of micro-organim swimming and other situations involving objects 
with moving boundaries. 

In the case of a spheroidal particle, the velocity disturbance v can be expressed as a perturbative 
series v = + ..., where is the value of v in the spherical case, given the same external 

flow conditions. It is standard practice to parametrize the deviation from spherical shape in terms 
of the quantity: AX (x) = X (x) — i?, where x = X (x) is the equation for the particle surface and 
the particle center is fixed imposing the condition that X does not have I = 1 components. The 
0{AX/R) contribution to v will then be in the form: 


= -AXd^v + = AXU' 


(3.7) 


plus contribution coming from possible 0{AX/R) corrections to U. Using Eqns. (2.9a-b) and 
(3.5), the term U' can be expressed in terms of the values U and U of the external shear and of 
the streaming velocity on the x = R surface: 


TT'e _ 21+1 ( SUtrn OTTe 'l i OlTTe 

Tjfm 2/+1 frm _i_ /+! frm 

^Im R ^Im R ^Im 


(3.8) 


From here, the 0{AX/R) contributions to force and torque can be expressed in a form which actu¬ 
ally does not require information on the detailed structure of the velocity disturbance; explicitating 
the 0{AX/R) corrections to U: 


F(i) = (F|UW) +^(r|AX|/iZm)C/' 

^Im 


Im 


and 


T« = (T|U(^’) +^(T|AX|/rZm)C/;^. 

^Im 


(3.9) 


(3.10) 


In fact, the matrix elements in Eqns. (3.9-10) require the calculation of integrals, which are basically 
equivalent to the ones needed in the reciprocal theorem approach |Q, say /^dA • cr • (U — U), 
where A is the particle surface defined by a; = X(x), and cr is the viscous stress generated around 
the particle when this is either translating or rotating in a quiescent fluid. 
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IV. Spherical particles with generic boundary conditions at the surface 


As a first application of the formalism illustrated in the previous section, the drag force and 
torque acting on a rigid, radius R spherical particle, carrying on respectively a translational and 
a rotational motion, can readily be calculated. In the case of a translation with velocity U, one 
finds from Eqn. (3.2): 

Fi = (F,|U) = -{3/2)fiRU j dn^ = -dn^iRU, (4.1) 

which is the standard Stokes formula for the drag. In the rotation case, take the angular velocity 
directed along 2 : 3 , with magnitude w. The velocity at the particle surface is then U = u}[—X 2 , Xi, 0] 
and the drag will be: 

Ta = (T 3 IU) = -3ujnR^ J dzd(j){l - z^) = -SttR^huj. (4.2) 

As a less obvious example, consider the velocity of a spherical particle at the centre of a quadratic 
shear flow: v = kx^gs. From Eqn. (3.2), the force is ^3 = — (FajU) = ^kR?2'k dzz^ = 2 ^ttkR^. 
Using the expression for Stokes drag, one obtains that the particle will move with a velocity 

kR? 

u = ( 67 r/ri?)“^F 3 = —^^ 3 . (4.3) 

The same result could have obtained using Faxen’s equation for the force |l7|| , and, in a certain 
sense, Eqn. (3.6) can be seen as an integral form of Faxen’s result, with the bonus that boundary 
conditions different from no-slip can be taken into account explicitly, and that perturbation theory 
becomes possible. 


A typical situation with moving boundaries is that of a swimming micro-organism. It is then 
possible to obtain an upper bound for the swimming efficiency, more precise than the one derived in 
1^. A definition for the swimming efficiency of a radius R spherical micro-organism is the following 
one, due to Lighthill Q: 


VL 


Qtt^Rv? 

W 


(4.4) 


where dirfiRu'^, from Eqn. (4.1), is the power necessary to set an identically shaped rigid particle 
in motion at speed u, and W = f dAf • U is the work per unit time actually carried on in 
swimming. From Eqns. (2.10a) and (3.5), in the case of a particle which does not change shape 
while swimming, this work is equal to: 


W = 


[{21 + 1 )+ 2i/L - UL)Ur:, + [{21 + 1)C7;- - (/ + 2)U;“ ) UZ*] (4.5) 


Im 


' \/K^ + 1 ) 







where U is the streaming velocity associated with the swimming stroke and U = —u. If u is along 
X 3 , it is easy to show that, necessarily: 


{FslelO) 

dnfiR 



-V^u!o 


(4.6) 


where the expression for the vector spherical harmonic Y®q = 1 ] — has been used 

[see Eqn. (2.2)]. In other words, only the jelO) component of U contributes to the force on the 
particle (there are no ”s” components, in any case, due to absence of particle changes of shape). 
From Eqn. (4.5) and using also the expression Y®q = obtains for the work performed 

in swimming: 


V2 

which sets the maximum swimming efficiency to 77 ^ = i, in the case the velocity associated with 
the ’’swimmimg stroke” is a steady jelO) circulation on the particle surface. This upper bound 
is stronger than the condition ryi < | presented in due to the neglect in Eqn. (11) of that 
reference of the contribution to dissipation from vorticity, which from Eqn. (2.8a) is given by: 


to-Uto)uro=l2^f^Rn^ ( 4 . 7 ) 



H = 






(4.8) 


As a last exercise, it is interesting to generalize the result provided by Eqn. (3.6) for the force on 
a spherical particle, to the case of an infinite surface tension fluid droplet. This clearly parallels the 
generalization of Faxen’s equation, obtained in using the reciprocal theorem. The boundary 
conditions for this problem are continuity of tangential stresses across the interface and absence of 
deformations from spherical shape, which means = 0. Focusing on the X 3 component of the 
force, it will be necessary to know only the jelO) component of U: 

F3 = -(E3|U) + (F3|elO)17ro (4.9) 

In order to calculate the velocity component U^q, the necessary boundary conditions are: 

U!o = U!o-Uto, Ul, = -Ulo and f!o = f!o + f!o (4.10) 

where /fg is given by Eqn. (2.10a) while /fg and /fg are given by Eqn. (2.10b). One obtains 
promptly: 
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where fi is the droplet fluid viscosity. Substituting Eqns. (4.11) and (3.2), into Eqn. (4.9) and 
using the expressions for Y^q previously given, one finds an expression for the force, again in the 
form provided by Eqn. (3.6): = — (Fi(/i//r)|U), where now: 

l^’*(A)) = + A)-i((A - l)ei + 5x,x) (4.12) 

Clearly, Eqn. (3.2) is recovered in the limit A = /t//r ^ oo. In the quadratic shear case, discussed 
in Eqn. (4.3), the force on the droplet will be: which is zero in the limit fi/fi ^ 0, 

typical of a spherical bubble. 

V. Spheroidal rigid particle in plane channel flow 
A. Geometrical approach 

The presence of solid surfaces bounding the flow leads to additional forces and torques on particles 
in suspension; these forces and torques can be obtained from Eqns. (3.6) and (3-9-10), once the 
modifications produced on the external velocity field v are known. If the distance I from the wall 
is large compared with the particle radius R, it will be possible to approximate the modification 
to V, with the first image at the wall, of the velocity perturbation by the particle. This image is 
obtained imposing on the wall the no-slip boundary condition v' -|- v = 0, where v is the velocity 
perturbation by the particle in the infinite domain case. 

If the particle is neutrally buoyant, the velocity perturbation at the wall will be v = 0{{R/l)‘^) 
and similarly, one will have that u" = 0{{R/l)'^) at the particle, while stress and vorticity will be 
0{{R/IY). In the far field regime, the main effect of the wall is therefore a force, which, as it is 
well known, will have a transverse component only if the particle is non-spherical. As illustrated 
in Fig. 1, a lot of what happens can be understood in terms of the particle and the flow symmetry 
properties, at least if the shear plane is also a plane of symmetry for the particle. 

Definite results can be obtained in the presence of one more symmetry plane, identifying a 
fore-aft line for the particle, or if the particle has fore-aft symmetry, if there are in total three 
symmetry planes. It is clear that in the two cases b and c, the lift will be proportional respectively 
to sin 20 and to cos0, where 0 is the angle of the particle symmetry axis with respect to the flow 
direction, as shown in figure. 

Transversal drifts arise also from the interaction of the particle non-spherical shape with the 
flow gradient. For instance, as shown in Fig. 2, a particle without fore-aft symmetry immersed in 
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Figure 1: Sketch of velocity perturbation around a particle in various shear flows, and image 
velocity field produced by a wall; the velocity fields are all in the particle’s reference frame, (a) 
perturbation around a spherical particle: the velocity lines are symmetric around the perpendicular 
to the wall (dashed line) and no lift is present; (b) perturbation by an ellipsoidal particle in linear 
shear: as in the spherical case, the field lines are characterized by a quadrupole symmetry; (c) the 
same for a particle without fore-aft symmetry in quadratic shear: in the near field, the velocity 
perturbation has an octupole symmetry. 

a linear shear flow will feel a force perpendicular to its symmetry axis. This means an additional 
drift proportional to cosO, analogous to the one described in Fig. Ic, in the case of a non fore-aft 
symmetric particle interacting with a wall, in a quadratic shear. In the perturbative case, this 
force arises from the combination in Eqn. (3.9) of the I = 2 harmonics in the external shear and 
the I = 3 components in AX. The same line of reasoning would lead to expect that an ellipsoidal 
particle in a quadratic shear should feel a force that becomes maximum when the particle axis 
is at 45° with respect to the flow. However, the issue is not so simple: given a flow in the form 
V = vo(a -|- b{x2 /the force will be in the form F = oF' -|- hY" with F' and F" in general 
not parallel. However, since the angle between the two force contributions is unknown, it is not 
clear what is the sign of Fj, when a is chosen in such a way to have Fa = 0. In fact, it will turn 
out in the perturbative regime, that F 2 = 0. 

Turning to the torque, similar symmetry based considerations allow to select which particle 
shapes and flow regimes lead to a non-zero contribution. In the perturbative regime, Eqn. (3.10) 
and the even parity of |F) require that AX and v be characterized both by even or odd spherical 
harmonics, which means, from the fact that \T) is / = 1, that these harmonics can only be equal. 
As with the force, this statement can be generalized to the strongly non-spherical regime to the 
result that quadratic shear can act with a torque only on particles whose shape is not fore-aft 
symmetric (see Eig. 3). Additional informations on the dependence of the torque on the particle 
orientation, can be obtained writing the external velocity field v in a reference system rotating 
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Figure 2: Forces on particles lacking fore-aft symmetry. In case (a), top velocity lines (going 
towards the right) do not correspond to bottom ones (going to the left); the end result is a force 
in the horizontal direction. To each line going upwards, instead, there is a symmetric one going 
downwards (see e.g. portions of the velocity lines to the extreme left and extreme right of the 
particle). Hence, no vertical component in the velocity is present. Analogous reasoning shows that 
in case (b), fore-aft asymmetry leads to a vertical force. 



Figure 3: Fore-aft symmetric particle in quadratic shear: no torques are produced. Because of 
symmetry of both the particle shape and the velocity field, for any surface element at r, there is 
another one at —r, such that the contributions to torque are equal and opposite. 
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with the particle. To fix the ideas, take v = v{x 2 )g^ and, indicating with a prime the rotating 
reference frame, take x'l = Xi and to be along the particle symmetry axis (hence, in Fig. Ic, 
the two reference frames coincide). One will have then in the case of a linear shear: 

03 X 2 = 63 X 2 cos^ 9 — 62 X 3 sin^ 9 — GgXg sin0cos6> + sin 0 cos 0 (5.1) 

while, for a quadratic shear: 

63 X 2 = 62 X 2 ^ sin 9 cos^ 9 + e^x'^ sin^ 9 — 203 X 3 X 3 cos^ 9 sin 9 

— 202 X 3 X 3 sin^ 9 cos 9 + 03 X 3 ^ cos^ 9 + 03 X 3 ^ cos 9 sin^ 9 (5.2) 

Consider first the case of a fore-aft symmetric particle in a linear shear flow. As illustrated in 
Fig. 4, because of symmetry across the planes x'^x^ and x^Xg the terms in sin 0 cos 0 in Eqn. (5.1) 
do not contribute. One recovers therefore, without having to calculate the velocity perturbation 



Figure 4: Derivation of the functional form for the torque on a fore-at symmetric particle in linear 
shear. The flow fields drawn in quadrants (a-d) correspond to the four terms, in the decomposition 
of the linear shear in the primed reference frame of 03 X 3 [see Eqn. (5.1)]. Elows (c) and (d) 
produce boundary conditions for the velocity perturbation v, which, because of their symmetry, 
do not lead to torque. 

around the particle, the expression for the torque, derived in [^ in the case of the ellipsoid: 

Ti =a(l-f 5cos26') (5.3) 

Turning to the case of a quadratic shear and of a particle which is not fore-aft symmetric, it appears 
from analysis of Eig. 5, that only the terms on the first line of Eqn. (5.2) contribute to torque, so 
that, in that case: 

Ti = c{l + dcos29) siii9 (5.4) 
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Thus, quadratic shear acts in the direction to align the particle axis with the flow. Conversely, if 
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Figure 5: Terms in the decomposition of a quadratic shear in terms of primed variables, which do 
not contribute to torque. Flows (a-c) correspond to the terms on the second line of Eqn. (5.2) 


linear shear is the only actor in the game and the particle is rigid, no equilibrium orientations are 
possible. (At least, this is true in the case of the rigid ellipsoid |^). It remains to be shown that 
quadratic shear plus fore-aft asymmetric shape, is actually the only mechanism for orientation 
fixing. In Fig. 6, it is shown how the last possibility of an odd contribution to torque from 
interaction between fore-aft asymmetry and linear shear is excluded. 



Figure 6: The torque on a fore-aft asymmetric particle in a linear shear, does not contain odd 9 
contributions. Case (b) is obtained from case (a) changing sign to the velocity at all points on the 
particle surface. Hence the torques in (a) and (b) are equal and opposite. Case (c) is obtained 
from (b) by a reflection. The velocity, and therefore also the force, remain the same, point by point 
on the particle surface. The position vector from xs to the point on the surface, however changes 
sign, and then the torque is the opposite to that in case (b). Hence the torque in (a) and in (c) 
are equal. 


To summarize, the istantaneous transversal drift of a particle in a channel flow is due to three 
distinct mechanisms, which, for shapes which have at least two planes of symmetry (one of which 
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is the shear plane), can be described as follows: 


i. An interaction with the walls due to the modifications in the velocity perturbation, produced 
by the particle shape deviation from spherical and the linear part of the external shear. This 
contribution is proportional to sin 20, with 9 the angle between the direction of the flow and 
the particle symmetry axis. This drift is expected to be directed away from the nearest wall, 
when the longer particle axis is aligned with the straining direction of the flow (see Fig. lb). 

ii. A second interaction with the walls, which is present when the particle shape is not fore-aft 
symmetric and the shear has a quadratic part. The resulting drift is proportional to cos0, 
and, at least in the near field regime (see Fig. Ic), one expects that it be directed away from 
the nearest wall, when the particle’s ’’tip” is in the direction of the flow. 

iii. A term proportional to cos 0, which is due only to the interaction of the particle with the 
gradient of v, and which is present only if the particle is not fore-aft symmetric. This drift 
pushes the particle towards the regions in which v is in the direction of the particle’s tip (see 
Fig. 2a). 


Using also Eqns. (5.3-4), this means that, in a channel flow, the law of motion for a particle away 
from the walls will be in the form: 


y = A{9, y, L) sin 20 -|- B{9, y, L) cos 9 
9 = a'(l -I- 5cos 20) -|- c'(l -I- dcos20) sin0 


(5.5) 


where y is the distance from a wall and 2L the channel width. From previous discussion, A and B 
are even under the exchanges 0 —> —9 and 0 —> 0 -|- m: 12, a' is proportional to the linear part of 
the shear, with |6| < 1; one has c ^ 0 only in the case of a fore-aft asymmetric particle in quadratic 
shear. From Eqn. (5.5) it appears that 0(0) = 0(7r — 0) and at the same time y{9) = —y{'K — 9). 
Thus, unless an equilibrium orientation 0o: 0(0o) = 0 exists for a' ~ c', y will have a periodic time 
dependence and no steady drift will occur. Of course, this result is a consequence of Stokes flow 


reversibilty |17 and is therefore valid also for generic non-symmetric particle shapes. This differs 
from the situation considered in where the deformability of the particle (a red cell) implied 
a contribution to the coefficient c' in Eqn. (5.5) proportional to cos0, meaning violation of the 
symmetry y{9) = —y{'K — 0) and the possibility of a mean drift also in the periodic regime (the 
region away from the channel axis corresponding to c' small). In the present case, instead, if the 
form of the coefficients A and B allows for a stable fixed point of Eqn. (5.5) at 0 = 0, y = 0, the 
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attraction basin of this point will be a smaller region close to the channel axis, and not the whole 
space domain. 


B. Analytical results 


The bra-ket approach of section III for the calculation of forces and torques can be generalized to 
the case of Vl, i.e. the part of the transverse drift due to the channel walls. From Eqn. (2.6a), for 
a neutrally buoyant particle, the far field behavior of v will be in the form: v = aijXiXj, 

with aij traceless and symmetric. Taking the channel walls to be perpendicular to the X 2 axis, 
analysis carried on in Appendix B shows that the particle drift in direction X 2 depends only on 0:22 
and can be approximated very well superposing the images produced separately by each wall. 
If the walls are located, with respect to the particle, at X 2 = Lr and X 2 = —Li, with Li,Lr > 0: 


0:22 / I 1 \ 

~\L^ ~ Zf/ 


(5.6) 


Following the approach of section III, one introduces a ket \L) such that: 022/4 = (T|v). Using 


Eqn. (2.9a) and the expressions for the spherical harmonics: F 20 ~ a\I f ~ 1) 42±2 = 

— X 2 ± 2 x 1 X 2 ), one finds: 022 = + \f^i ^22 + ^ 2 ,- 2 )^ terms involving 

scalar components which will be shown later not to contribute. From Eqn. (2.4) and the 
expression given above for y 2 m) one finds therefore: 


\L) = (1 - xl)x 2 , -xlxs] (5.7) 

plus the contributions associated with C//^. From here, the generalization of Eqns. (3.6) and (3.9) 
can be obtained, but the first one gives identically zero, and we are left with: 


^ ^ film 

and the absence of /i = s contributions in the sum [see Eqn. (3.8)] guarantees that the 
contributions could actually be neglected in the calculation of |L). 

The law of motion of a particle in a channel flow can then be determined perturbatively using 
the formalism provided by Eqns. (3.2), (3.4), (3.8-10) and (5.7-8). Indicate the transverse particle 
coordinate in the laboratory frame with y, and fix the origin of the axis at the left wall, so that 
y = Li; in this reference frame, the velocity profile is ujy{l — y/(2L)) where 2L = Li + Lr- As 
before, {xi,X 2 ,X 3 } is the reference frame translating with the particle center, with X 3 along the 
flow and X2 along the direction of variation for v; for a particle at y, from Eqn. (4.3): 


V = (w(l - y/L)x 2 + {uj/L){R^/3 - xl))e 3 . 


(5.9) 
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To fix the ideas, consider an axisymmetric particle whose shape in the co-rotating frame {x'l , a;^} 

= {xi , X 2 cos 9 + x^, sin 0, x^, cos 6 — X 2 sin 9} is: 


AX(x)/i? = e(-P 2 (i^) + aPsix’s)) 


(5.10) 


with e small, e,a > 0, and P 2 {z) = (3/2)z^ — 1/2 and Ps^z) = {5/2)z^ — (3/2)z Legendre polyno¬ 
mials. Thus, for 0 = 0 the particle is an oblate ellipsoid, and, for a > 0, the tip of the particle is 
at a ;3 > 0. (By composition of the harmonics in v and in AX, only I = 2 and I = 3 terms in AX 
give non-zero contribution to the matrix elements for force, torque and wall produced lift). 

The particle instantaneous transversal velocity is the sum of the contribution from the wall, 
given by Eqn. (5.8), and the one due to interaction between the fore-aft asymmetry of the particle 
shape and the shear: y = Vl+ Substituting Eqns. (5.9-10) into Eqn. (5.8), and using Eqns. 

(3.8) and (5.6) gives Ul- Repeating the same operation with F 2 , with Eqns. (3.2) and (3.9) in 
place of (5.7-8), and using Eqn. (4.1): 


y = eujR 


(1 - y/L) cos^ ^ “ 7 ) ^ 


/ 3 , , , ^ 25aR /5 , ^ \ A f 

+(55(1 - »/i)».2e + —(jcos e- 1) coso) A 


R^ 


(5.11) 


■r i^L-y)^. 

Details of the calculation are given in Appendix A; see in particular Eqns. (A14-15). Notice that 
no contribution from the quadratic part of the shear is present in the term in the first line of the 
equation: the one which doe not come from interaction with the wall. Anyway, this seems to be 
fortuitous and non-zero contributions may show up at higher orders in the perturbative series in e. 

The contribution to torque from AX is obtained working with Eqns. (3.4) and (3.10) and 
this leads to the equation of motion for the particle orientation, using also Eqn. (4.2): 9 = 
1(1 - y/L) + 87 ^- Hence, from Eqn. (A16): 


/I 

3e „ A 

R, 

H5 2/, 

9 \ 

( + 

— cos 29 ) 

— eau !— 

— cos^ 9 — 

TT 

V2 

4 / 

L ' 

y 8 

14/ 


(5.12) 


Thus, close to the axis (i.e. when \L — y\ < eaR), a particle with fore-aft asymmetric shape will 
be forced to get the orientation 0 = 0. Under these conditions, the dominant part of the drift is 
the interaction with the shear [first line of Eqn. (5.11)], which pushes the particle towards y = L. 
Away from it, the standard Jeffery solution in the small eccentricity limit corresponding to a 
periodic motion regime, is recovered, and, as discussed at the end of the previous sub-section, the 
mean transverse drift is equal to zero. 
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VI. Conclusions 


The present work has been carried on along two parallel lines: the derivation of qualitative 
results on particle behaviors, by means of simple geometrical reasoning, and actual quantitative 
calculations, using perturbation theory around the spherical particle case. Emphasis was put, 
especially in the geometrical proofs, on the origin of the results in the symmetries of Stokes dy¬ 
namics. It was thus possible to obtain a classihcation of the behaviors that a particle, depending 
on its shape, can have in linear and quadratic bounded shear flows. The present analysis was 
limited to shapes with at least two planes of symmetry, so that a large number of more irregular, 
strongly non-spherical objects were excluded. Also, some of the results, namely Eqn. (5.5) rested 
on a stability hypothesis for configurations, in which one of the particles symmetry planes coincide 
with the shear plane. (It should be mentioned that recent results suggest instability of such 
configurations, and also chaos, at least in the case of very elongated ellipsoids). In spite of these 
limitations, it was possible, in several important cases, to give a more direct physical content to 
results about the functional form of torques and forces, without actually having to solve Stokes 
equation; one example is the generalization of Jeffery’s result for the torque on an ellipsoid JT^ , 
provided by Fig. 4. 

One of the motivation for the present study was to get some ideas on the mechanisms deter¬ 
mining transverse migration in plane shear flows. Some of the results obtained, were implicit in 
previous works on the behavior of droplets [l^ and vesicles |p^ , in viscous shear flow. However, 
an explanation was missing, of which ingredients in the dynamics of these systems, were at the 
basis of the migration process. One result has been to clearly illustrate the role of quadratic shear 
and of particle deformability. Summarizing, a rigid particle will drift to the center of a channel 
flow, only if fore-aft asymmetric and sufficiently close to the channel axis, for quadratic shear to 
be dominant. A particle will drift to the channel axis irrespective of its initial position, only if 
deformable. This may happen, either because of deformations, produced by the quadratic shear, 
which break fore-aft symmetry in a preferential direction |^, or because of the presence of tank- 
treading motions |^, |^), or because of a combination of the two effects, as in the case of droplets 

0 - 

The goal of this paper was, from one side, to understand the origin of particle behaviors in their 
shape symmetries and in those of the fluid flow; from the other, it was to furnish an efficient method 
for carrying on calculations in the perturbative quasi-spherical regime. Apart of the illustrative 
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examples considered in section IV, the technique allowed to obtain quantitative informations on the 
transverse drift of axi-symmetrical particles in channel flows. As clear already from the geometrical 
considerations presented before, a particle drifts basically because of three mechanisms. These are 
the interaction of the I = 2 harmonics in the particle shape with the wall, mediated by the linear 
part of the shear; the direct interaction of the I = 3 shape harmonics with the linear part of the 
shear (even in the absence of walls); the interaction of the I = 3 shape harmonics with the walls, 
mediated by the quadratic part of the shear. For a rigid particle in a channel flow, only the second 
contribution plays a role, given the absence of steady drifts close to the walls. For different reasons, 
the same situation is met in the case of quasi-spherical vesicles in channel flows . 

One purpose of the bra-ket formalism described in Eqns. (3.9-10) and (5.8) was to provide 
a viable tool in situations where a reciprocal theorem approach is too cumbersome. In the solid 
particle cases, considered in this paper, the amount of work needed in the analysis, was clearly 
equivalent to what would have been necessary in a reciprocal theorem approach. Thus, the choice 
of the technique was mainly a matter of taste. The situation becomes different if the fluid is 
strongly non-Newtonian, or in general, if one has to deal with constitutive laws of an unusual form. 
In a separate paper, the present technique was used to study the dynamics of inextensible 
membranes. In that case, the presence of constitutive laws different from incompressibility made 
an approach based on ordinary perturbation theory, rather than on the reciprocal theorem, the 
only reasonable choice. 

Aknowledgements: I would like to thank Howard Stone and Dominique Barthes-Biesel for inter¬ 
esting and helpful conversation. Part of this research was carried on at CRS4 and at the Laboratoire 
de Modelisation en Mecanique in Jussieu. I would like to thank Gianluigi Zanetti and Stephane 
Zaleski for hospitality. 

Appendix A. Calculation of vector spherical harmonics components 

Matrix elements involving vector spherical harmonics, in the form given by Eqn. (2.3) can 
be expressed in terms of integrals containing only plain spherical harmonics. For instance, it is 
possible to write: 

f dB,r(a:Vy;);,(x*)) • v(i?x) 

+ 1))2 Jo J 

= {l{l + 1 ))“^ J dBa;Y)^(ii)[3x • v(i?x) - V • (a:v(i?x))]x=_R (Al) 
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which is the second of Eqn. (2.4); one obtains the third equation in that system in completely 
analogous manner. For ease of notation, it is possible to introduce operators such that Eqn. 
(2.3) read: 

uL(i?) = (H 6 L-v) (A2) 


where (Imlf) = f dflxYi^{xi) f (x) indicates a plain spherical harmonics component. 

For the analysis in section V, it is necessary to calculate quantities in the form • v with 

V = [0,0, 0 : 2 ], [0,0,1/3 — X 2 ] and 6^^ ■ AX|Z), where Z = F 2 ,Ti,L and AX is given by Eqn. 
(5.10). Taking for simplicity i? = 1 from now on, one obtains promptly, in the linear shear case: 


OL • V = 


3X2X3 

\/i(r+T) 


Xi 


\/K^ + 1 ) 


while, in the quadratic case: 

6 


V = 


2 x 3 (xf — 5 x 2 + ^ 3 ) 


2xiX2 


(A3) 


(A4) 


d^JUTTY) ™ ^/W+Y) 

From here, the velocity components for the unperturbed flow are obtained. For the linear shear 
case: 

= ii/TT/lS, = ii/VS and = ±^/YJz (A5) 

In the quadratic shear case: 

4C^ 2C^ 2C^ 4C2E/ 

•^3 

2y/TYK 


TP — 
*^10 — 


45 


TP — — 

" 45 


TT^ — 
^30 — 


TP — 

'^3,±2 — 


27r 

105’ 


35 


TP — 

*^30 — 


105 


Ul±2 = 


105 


- 2iV^ 

and C^ 2,±2 “ T 


(A 6 ) 


The calculation of the terms 6 ))^ • AX\Z), necessary for the matrix elements {Z\AX\ij],m) entering 
Eqns. (3.9-10) and (5.8) is more cumbersome. Life becomes simpler, however, exploiting the 
knowledge of the dependence on the orientation 9 that the resulting forces, torques and drifts 
must have. To make an example, as discussed in section V, the contribution from P 2 {z') to Ul, 
is necessarily proportional to sin26l. Isolating that contribution to P 2 {z'), i.e. —( 3 / 2 )x 2 X 3 , one 
obtains for instance: 


0 


Im 


P2\L) 


45 X 2 X 3 ( 2 xf — 3x| -|- 2 x 3 ) 


(A7) 


sin 29 647r y/Kl + 1 ) 

plus terms that will give zero after integration in dfla,. The hnal expressions for the non-zero 
matrix elements {Z\Pa\^lm) are given below. For the lift from the wall: 


(L|P 2 |e 2 , ± 1 ) = sin 2 d, (L|P 2 |ml ± 1 ) = ± A sin 2 d 


(A 8 ) 
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and 


(LiPslelO) = 1 X 2 y^(5 cos^ 0 — 4) cos0, (LlPajeSO) = y^^(5 cos^ 0 — 4) cos0, 

(L|P3|e3,±2) = 4igY^(5cos2 0-3)cos0, (L|P3|m2, ±2) = y^( 4 - 5 cos^ 0) cos0. 

For the force, indicating |/) = (67r^i?)“^IF 2 ): 

(/|P3|e2,±l) = ^y^(4- 5cos2 0)cos 0, (AlO) 


(/|P 2 |elO) = -^Y^sin20, (/|P2|e30) = -^Y^sin20, 

(/|P2|e3,±2) = (/|P2|m2,±2) = ±^^sin20. 

Finally for the torque, indicating \t) = (87r/rP^)“^|Ti): 

(t|P 2 |e 2 ,±l) =T^tpcos20, (t|P 2 |ml,±l) = 

40 V TT 40 V TT 


(All) 


(A12) 


and 

3 3 fro 

(tiPsleSO) = —^\ —(5cos^0- l)sin0, (t|P3|e3,±2) = —\ —(4cos^0- l)sin0 
iiz V "TT oyo V TT 

(t|P3|m2,±2) = =F-^i/^sin0. (A13) 

i4U V TT 

The expressions for lift, force and torque are obtained substituting Eqn. (A5-6) into Eqn. (3.8) 
and then, together with Eqns. (A8-13) into Eqns. (3.9-10) and (5.8). If no torques act on the 
particle, the particle will spin in such a way that in the third of Eqn. (3.8). 

Using Eqns. (5.8) and (A8-9), one obtains for the wall contribution to lift, respectively, in the 
two cases AX = P 2 {z'), v = [0,0, X 2 ], and AX = Psiz'), v = [0, 0,1/3 — x|]: 

Vl =-— f{Li,Lr) and Vl = — (^-cos -10j cos0/(P/, Pr-) (A14) 

where f{x,y) = R^{x~^ — y~‘^) [see Eqns. (5.1) or (B26)]. Next, from Eqns. (3.9) and (AlO-ll), 
one finds for the drift produced by the force P 2 (i.e. the contribution to drift that is independent 

from the wall), in the two cases AX = P 2 iz'), v = [0,0,1/3 — x|] and AX = Psiz'), v = [0, 0, X 2 ]- 

Vf=0 and u/= — y cos^ 0^ cos0 (A15) 

where Vf = (67r/iP)“^P2. Einally, from Eqns. (3.10) and (A12-13), one finds the 0{AX) correction 
to the particle angular velocity 0^^/ in the two cases AX = P 2 {z'), v = [0, 0, X 2 ] and AX = Psiz'), 
v= [0,0, l/3-a;|]: 

0^"^' = --cos20 and 0*^' = ( -cos^ 0^ sin0 (A16) 

4 Vl4 8 / ^ ^ 
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Appendix B. Wall contributions to transverse drift in channel flows 


Consider a fluid flow, bounded by two walls perpendicular to the X 2 axis, located respectively 
at X 2 = Lr and X 2 = —Li = — 2L. We want to calculate the lift along the X 2 direction, on a 

small particle placed at the origin of the axes, in the regime i? <C Li, Lr, where R is the size of the 
particle. In this regime, if the the flow lines are taken parallel to the walls, the lift will be given 
by the image at the walls of the far field part of the velocity perturbation. 

It is therefore necessary to know the velocity in the interval —Li < X 2 < Lr, given boundary 
conditions at X 2 = Lr, X 2 = —Li. It is expedient to Fourier transform in xi_ 3 . The vorticity and 
continuity equations, given stationary, creeping flow conditions lead to the system: 

{ 82 X 3 — ik3V2 = aexp{k{x 2 — Lr)) + bexp{—k{x 2 + Lij) 
iiksvi - kivs) = f exp{k{x 2 - Lr)) + g exp(-fc(a; + L/)) 
ikiV 2 — 82 X 1 = cexp(A:(a ;2 — Lr)) + fiexp(—A:(a: + Li)) 

82X2 + i{kixi + /cat’s) = 0 . 

where k = y^kf + fc|. Multiplying the first and third equations in this system by fcs and ki and 
taking the difference, one gets: 

82 i{kiXi + ksXs) + k‘^X 2 = /(fcso — kic) exp(fc(a ;2 — Lr)) + 1(^36 — kid) exp{—k{x 2 + Li)) (B2) 


and using the continuity equation: 

(5| — fc^)v 2 = —/(fcao — fcic) exp(fc(a ;2 — Lr)) — iiksb — kid) exp{—k{x 2 + Li)) (B3) 
which leads to the expression for X2' 

X 2 = (^X 2 + iA{Lr — X 2 )'j exp(k(x 2 — Lr)) + (^X 2 + iBiLi + X 2 )^ exp{—k{x 2 + Li)) (B4) 
where: 

ksa-kic ksb-kid 

A = and B = (B5) 

Substituting into the first and third of Eqn. (Bl) and integrating in X 2 ' 

X 3 =V3 + + ^(a;2 - Lr)"j exp{k(x2 - Lr)) 

- (^^4^ - W - ^(ac2 + ii)) expi-k{x2 + Li)) 

Xi =■01 + W + ^(^2 - Lr)j exp(fc(x2 - Lr)) 

- - hiB(^x 2 + Li)'j exp{-k{x 2 + Li)) 

The integration constants Oi^a are eliminated using the second and fourth of Eqn. (Bl). Noticing 
the fact that the left hand side of those equations must be combinations of terms proportional to 
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exp(±A:a: 2 ), we get: 


k 3 ,vi — kiV 3 = 0 
kiVi + ksvs = 0 


(B7) 


leading to: t)i ,3 = 0. 

We now can enforce boundary conditions on Vi^^: 


Uo = 


_ / ik^V 2 -\-a 


( 


ksA 

k ~W 

ik3V2~^b _ k^B 
k k'^ 

—C 


.,r _ (ikiv^-c kiA\ , f ikivl^-d , 

^ + 


) + ( - + ^(1 + 2fcL)) eM-^kL) 

) + ( - ^^4^ - ^(1 + ^kL)) exp(-2fcL) 


k k'^ 

2kL)) exp(-2A:L) 

ikivL—b 


(B8) 


^ ^ I' _ ^k^^ - ^(1 + 2fcL) j exp(-2fcL) 

Adding the right and left terms, with factors ki and ks respectively for components 1 and 3, we 
get: 


I A + { 2 kL — 1 ) exp{— 2 kL)B = —ikv2 + ikv2 exp(— 2 fcL) + k^v^ + kivl 
+ { 2 kL — 1) exp(— 2 A:L)A = —ikv2 + ikv2 exp(— 2 fcL) — k^v]^ — kiv\ 

Finally, the boundary conditions on V2 are obtained from Eqn. (B 4 ): 

jv2=V2 + {^2 + exp(— 2 A:L) 

1 ^2 = ^2 + (^2 + ‘ 2 -iLA) exp(— 2 fcL) 


(B9) 


(BIO) 


At this point, the coefficients A, B and can be obtained from Eqns. (B8-9): 


v^2 = c(8ikBe-^’^^{kivl + k3v;) + il-{l-4kL + 8 {kLf)e-^’^^)vl 2 

+2iLe-'^^^{l - e-^'^^){kiv[ + fcgu') - (1 + 2kL + {2kL - (Bll) 

vl = c(^-8ikBe-'^'^^{kiv[+k3V^3) + {l-{l-^kL + 8{kLf)e-‘^'^^)vl 

-2iLe-2'=^(l - e-'^'^^)ikivl + k3vl) - (1 + 2kL + {2kL - (B12) 


A = c(^(l-(l + 4kL)e-^'^^)(kiv^ + k3V^)+i((l-4kL)e-^^^-l)kv^ 

+ ((4kL - 1 + e-^^^)e-^^^(kivi + + t(l + 4fcL - (B13) 

B = c{^-(l-(l + 4kL)e-^’^^){kiv{ + k3v[)+i{{l-4kL)e-‘^’^^-l)kv ^2 

-{{4kL - 1 + {kiv{ + /cguj) + i(l + 4fcL - (B14) 

where 

C = (1 - (2 + 16(fcL)2)e-‘‘'=^ + (B15) 

Erom Eqn. (B4), to calculate V 2 and therefore the lift, we do not need the coefficients a — d; 
knowing A, B and is sufficient to the goal. Using Eqns. (Bll-15), we obtain: 


V 2 = iF{Lr, k){kivl + k3vl) + G{Lr, k)v 2 — iF{Li, k){kiv[ + k3v\) + G{Li, k)v 2 (B16) 
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where: 


and 


F{Ly, k) = C{Ly - {Ly + SfcL^ _ 4/cLLy)6“”2^) 
-{Ly - 8kL^ + AkLLy)e-^'^^ + 


(B17) 


G{Ly,k) =C{l + kLy-{l + k{AL-Ly)+AeL{2L-Ly))e-^<^^^-^^y^ , . 

-(1 - k{AL - Ly) + Ak^L{2L - Ly))e-^^^ + (1 - kLy)e->=^^^-‘^^y'))e-'^^y 


The lift is obtained evaluating the inverse Fourier transform of V 2 with respect to xi and 0 : 3 , at 
X = 0: 


^ pOO pOO p 27 T p 27 T 

Vi, = —-,—TT / kdk / xdx / dtp / di^expl—fcicos^/>) 

( 27 r)^ Jo Jo Jo Jo 

X ^ (±iF(Lj^, fc)(A;iwf (x) +/c 3 w|(x)) + G(Lj^, A:)-u|(x)) (B19) 

y=l,r 

where x = xi = x cos{tp + cp) and ki = k cos cp. The ± signs go respectively with y = I 

and y = r. As boundary conditions on the wall, we have v'^ = — v, where v is the far field part of 
velocity perturbation by the particle. In the absence of external forces, this is in the form: 


v(x) = X ^((a — 0)xl + pixl — axl)x 


(B20) 


plus terms which are not symmetric in Xi and X 3 , and give zero after integration. Substituting 
into Eqn. (B19) and carrying out the integrals in (p and tp, we get, in terms of Bessel functions 
Jy{x): 


g fOD ^00 2^2 _ ^2 

Ul = - / kdk xdx ±{G{Ly,k)LyJo{kx) - F{Ly,k)kxJi{kx)) ^ r9.\3/2 (^^1) 

2 Jo Jo 

where the following formulae for Bessel functions have been used ||2^ : 

p 2 'K p 2 'k 

/ d^/>exp(—iacosi/;) = 27rJo(a) and / di/^cosi/;exp(—iacos^/>) = 27riJi(a) (B22) 

Jo Jo 


P3 1 1 


(B23) 


The i-integration can be carried out analytically as well, using the formulae 


xJo{kx)dx 1 F kL 
lo (x2 + L2)f “ 3L3 


-kL 


^ ^3 


Jo{kx)dx 2 — kL 


-kL 


0 (X2+L2)§ 


3L 


roo 2 


x^Ji{kx)dx ke 


-kL 


lo (x2+L2)l 


3L 


and 


a;^Ji(fca:)dx , kL. 
lo (a:2 + L2)| 


(B24) 


24 












We thus get the hnal result: 



fc^dfc Y, MFiLy, k){l - kLy) + G{Ly, k)Ly)e-'^^y 
V=l,r 


(B25) 


and, for oo, using the asymptotic expressions F{Li,k) = Lie and G{Li,k) = (1 + 

kLi)e~^^’-, this leads to the expression for the lift obtained in [p0|: 


JL 

4Lf 


(B26) 


[One factor 1/2 was missing from |^, Eqn. (51), due to a calculation mistake there]. In the 
general case, the remaining integral in k must be carried out numerically. In Fig. 7 below, the 
dependence of Ul resulting from Eqn. (24) is compared with the analogous expression obtained 
using the large Ly asymptotic expressions for F and G: 



Figure 7: Behavior of the lift Ul in function of the distance Li of the particle from the channel 
left wall. The results of Eqn. (B24) and of its approximated version Eqn. (B26) are almost 
undistinguishable 
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